Theoretical Results for Sandpile Models of SOC with Multiple Topplings 
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We study a directed stochastic sandpile model of Self-Organized Criticality, which exhibits recur- 
rent, multiple topplings, putting it in a separate universality class from the exactly solved model of 
Dhar and Ramaswamy. We show that in the steady-state all stable states are equally likely. Then 
we explicitly derive a discrete dynamical equation for avalanches on the lattice. By coarse-graining 
we arrive at a continuous Langevin equation for the propagation of avalanches and calculate all the 
critical exponents characterizing them. The avalanche equation is similar to the Edwards- Wilkinson 
equation, but with a noise amplitude that is a threshold function of the local avalanche activity, or 
interface height, leading to a stable absorbing state when the avalanche dies. It represents a new 
type of absorbing state phase transition. 
PACS numbers: 05.65.+b, 87.23. Ge, 87.23.Kg 



I. INTRODUCTION 

Sandpile models of stick-slip dynamics have received 
considerable attention as canonical models of self- 
organized criticality (SOC) pj. SOC refers to the 
widespread tendency of many extended, dissipative dy- 
namical systems to evolve inevitably towards a complex 
state with power-law correlations in space and time: a 
"critical" state. Of course, a critical state is only one pos- 
sible example of complex phenomena that can emerge in 
large, self-organizing systems composed of many strongly 
interacting parts. No doubt there are other types of com- 
plex states that have not yet been so well characterized 
mathematically, e.g. for example in networks S. From 
this viewpoint, the phenomena of SOC itself is a proto- 
type for how complexity emerges in nature without fine 
tuning parameters. In spite of the gross simplicity of 
various cellular models that have been introduced, and 
hundreds if not thousands of numerical studies of SOC, 
only minimal analytic understanding has been achieved. 

In fact, a survey of analytic works on sandpile mod- 
els of SOC is exceedingly short. The model of SOC in- 
troduced by Bak, Tang, and Weisenfeld (BTW) f| has 
yielded to some analytic treatment associated with its 
abelian properties, primarily due to the work of Dhar and 
collaborators 0]. The scaling properties of waves, where 
each site only topples, or releases grains, once has been 
understood by Priezzhev and collaborators ^|J|]. Nev- 
ertheless, the large scale properties of avalanches, where 
each site can topple many times in response to a single 
grain being added to the system, remain unsolved and 
the numerical situation controver sial [|-|. The same is 
true for the Zhang model where some limited progress 
has been made using methods from dynamical systems 
theory E|. Recurrent, multiple topplings within an 
avalanche also appear in most other unsolved sandpile 
models, such as the stochastic Manna model [[nj, the 
universality class O] represented by the Oslo rice pile 



model |12|, cellular models of vortex dynamics jl4| ; as 
well as one-dimensional trough models exhibiting multi- 
scaling jl5| . The difficulties preventing progress in solv- 
ing any of these simplified models in particular, or find- 
ing general analytic tools for granular systems exhibiting 
SOC, appear to be related, in part, to the existence of 
recurrent topplings. 

This statement is further supported by the following 
facts: Dhar and Ramaswamy (DR) |ll| introduced a di- 
rected version of the BTW model, and solved for the 
avalanche distribution and many other properties ex- 
actly. In the DR model, it can be rigorously proven that 
no multiple topplings occur. (Consequently, the elegant 
DR solution, as it has been conceived thus far, does not 
address the full complexity of discrete or granular mod- 
els of SOC.) The fixed scale transformation method of 
Pietronero and collaborators |l7j also explicitly ignores 
the presence of multiple topplings. One consequence of 
this fact is that this method puts the stochastic Manna 
model and the BTW model into the same universality 
class, which is not consistent with most numerical works 
|t],[L8| (except || ) , including those measuring unequivocal 
differences in aging behaviors Multiple topplings, 

where the activity can return an arbitrarily large num- 
ber of times, do not appear in any mean field description 
p0| , since in high enough dimensions, the avalanche ac- 
tivity is not recurrent, or able to return more than a finite 
number of times at any site. Multiple topplings are 
a fluctuation effect associated with self-intersections of 
the avalanche cluster in space and time Q , and as such 
are relevant below some upper critical dimension. 

Certainly, the intricacies associated with multiple top- 
plings are not the only ones that present themselves in 
attempting an analytic treatment of granular models of 
SOC. For example, the fact that the dissipation process 
is confined to the boundary, which forces the system to 
self-organize, is an important and subtle point because 
the boundary cannot be scaled out in the limit of large 
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system sizes as is usually done in statistical physics. In 
principle, the boundary is always important, because the 
incoming sand grains must be transported to it, no mat- 
ter how large the system size. The broken translational 
invariance associated with the boundary often leads to 
long range boundary effects in the metastable states (see, 
for example, [^3|,|^|). It might be useful to pry these com- 
plications apart, treating one issue at a time. Here we 
focus on the problem of recurrent or multiple topplings, 
and seek a model which does not present other difficul- 
ties. 

Recently, Pastor-Satorras and Vespignani p5j have 
studied numerically a stochastic directed sandpile model 
(SDM), which is a stochastic version of the exactly solv- 
able model introduced by DR. This stochastic model is 
simpler and presumably unrelated to the directed mod- 
els introduced and studied by Tadic and collaborators 
|p6| . Pastor-Satorras and Vespignani demonstrated by 
numerical simulations that the model exhibits multiple 
topplings which changes the universality class, making it 
distinct from the DR model. This was accomplished by 
numerically measuring and comparing various critical ex- 
ponents characterizing the avalanches. Its close relation 
to the DR model, which has an exact solution, suggests 
to us an analytic study. 



as the temperature. The height of the interface represents 
the number of topplings in an avalanche. The steady 
state that is eventually reached in the limit of large times 
is always the state of no activity where the height of the 
interface is zero everywhere and the avalanche has died. 
Thus the equation describing avalanche dynamics corre- 
sponds to an absorbing state phase transition where the 
the transient state is governed by the EW equation in 
the region where it survives. Section V also describes 
an analysis of this nonlinear equation. We extract all 
the (nontrivial) critical exponents for avalanches, i.e. in 
d = 1, D = 7/4, r = 10/7, z = 2, r t = D = 7/4 distinct 
from the DR model. For d > 2, where multiple topplings 
are not relevant, the critical exponents are the same as in 
the DR model. All of these results agree perfectly with 
previous numerical works [§5],^8) . We also write down the 
Fokker-Planck equation for the probability distribution 
of the number of topplings at each site in an avalanche, 
although we do not solve it. Finally, we conclude with a 
brief comment on some possibilities for future analytical 
work on absorbing state phase transitions and granular 
models of SOC. 



II. DEFINITION OF DIRECTED MODELS 



A. Summary 

We proceed with an analysis of the SDM as follows: 
First we define the DR model and the SDM. For peda- 
gogical reasons, in Section III, we review the proof that 
the critical state of the DR model is the set of all sta- 
ble states with equal probability. We also review some 
necessary parts of Dhar's construction of an operator al- 
gebra for stochastic models, which can be represented as 
deterministic models with a quenched array of random 
numbers. Combining these two works, we then show that 
for the SDM the critical state is also the set of all sta- 
ble states with equal probability, described by a product 
measure. Using this fact, we show in Section IV that the 
SDM can be recast as a generalized branching process 
propagating in an uncorrelated environment, enabling a 
study of the infinite system. By carefully analyzing the 
microscopic dynamics of this process on the lattice, we 
explicitly derive a discrete dynamical equation for the 
propagation of flowing grains in avalanches. In Section V, 
coarse graining this discrete equation gives a continuum 
equation for avalanches that should describe the large 
scale properties of any microscopic model with the same 
symmetry, conservation of grains, and stochastic effects. 

Notably, our equation is similar to the Edwards- 
Wilkinson (EW) equation except that the amplitude 
of the nonconservative noise is a Heaviside (theta) func- 
tion of the local activity. Crucially, the noise amplitude 
is a threshold function, rather than being a constant, such 



Consider a two dimensional square lattice as shown in 
Fig. 1. The direction of propagation is labeled by t, 
with < t < T. The transverse direction is labeled by 
x, with periodic boundary conditions. Only sites with 
(x + t) even are on the lattice, so that a; is a positive in- 
teger modulo 2X, and the lattice has a total of TX sites. 
On each site, an integer variable zix, t) is assigned. The 
i'th grain is added to a randomly chosen site Xi on the 
top row t — 0. There z(xi,0) —* z(xi,0) + 1. When any 
site acquires a height greater than z c — 1 it topples, i.e. 
z(x, t) — ► z(x, t) — 2 for z(x, t) > z c . 
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FIG. 1. Directed Sandpile Models Grains from active 
sites in row t topple onto sites in row t + 1. For example, 
grains from the site indicated by the open circle topple only 
onto the nearest neighbor sites indicated by the arrows. 

The two models differ with respect to the transmis- 
sion of grains out of a toppling site. In the DR model, 
one grain is transferred to the left downstream neigh- 
bor and one grain to the right so the toppling rule is for 

z(x, t) > z c 



z(x, i) 

z(x- l,t + l) 
z(x+ l,t + l) 



z(x,t) - 2 
z(x-l,t + l) + l 
z(x+l,t + l) + l 



For the SDM, on the other hand, each grain from a top- 
pling site is given equal probability to go to any down- 
stream nearest neighbor. In this case, when the site (ar, t) 
topples, 



z(x, t) —* z(x, t) 



and 



z(x - l,t+ 1) - 
z(x+l,t+l) - 

with probability 1/2, or 

z(x -l,t + l) 

z(x + l,t + 1) 

with probability 1/4, or 

z(x- l,t+ 1) - 
z(x+ l,t + 1) - 



z(x-l,t+l) + l 
z(x+l,t+l) + l 



z(x-l,t+l)+2 
z(x + l,t+ 1) 



z(x- l,t + l) 
z{x + l,t + l)+2 . 



with probability 1/4. Thus, the SDM is a directed ver- 
sion of the model introduced by Manna. 

In both directed models, grains are conserved during 
each toppling event. This is true except at the open 
boundary t = T where toppling sites simply discharge 
their grains out of the system. Sites are relaxed accord- 
ing to a parallel update until there are no more unstable 
sites, and the properties of the resulting avalanche are 
recorded. Then a new avalanche is initiated by adding 
a single grain to a randomly chosen site on the top row, 
t = 0. An avalanche can be characterized by its longitu- 
dinal extent, t c , the largest t row affected, its width, x c , 
the largest transverse distance from the avalanche origin 
to any site affected by the avalanche, its area, a, the total 
number of sites affected, its size, s, the total number of 
toppling events, and the maximum number of topplings 
at a site, n c . 

It is straightforward to generalize this definition to 
higher dimensions, with the number of directions trans- 
verse to the direction of propagation being d. In this 



case z c = 2d — 1. At a toppling site z — > z — z c — \. 
In the DR case each downstream neighbor receives ex- 
actly one grain. In the stochastic case, each downstream 
neighbor has equal probability l/2d to receive each grain. 
For simplicity of notation and concepts we will focus our 
discussion on the case d = 1 unless otherwise noted. 



III. STATES ON THE ATTRACTOR 

For both directed models, any configuration satisfying 
< z(x,t) < z c for all (x,t) is stable. The total number 
of such configurations is (z c + 1) TX . For clarity, we now 
review the argument showing that in the steady state, all 
such stable states are equally likely in the DR model. 



A. Review of Some Exact Results by Dhar and 
Ramaswamy 

Let Co be a starting configuration with the i'th parti- 
cle added at site Xi , resulting in the new stable configura- 
tion Ci . Then Cj is uniquely determined by the dynamics 
given Xi and Cj_i. The crucial point is that this dynam- 
ics is invertible. On the top row Cj differs from Cj_i 
only at the site Xi, with z(xi,0) in Cj being more than 
its value in C_i by one (mod2). Other rows in C_i 
are the same as in Ci if there was no toppling at (xi,0); 
otherwise the z's in the first row, t = 1, in C,_i are the 
same as in Ci, except at the two downstream neighbors, 
{xi — 1, 1) and {xi + 1,1) of (a;*, 0) whose heights are less 
by one (mod2) than their values in C,-. This obviously 
continues for subsequent rows. Thus given Ci and Xi we 
can uniquely determine C_i. 

For a given Ci, there are precisely X distinct choices 
of Ci-i and C+i corresponding to X distinct, possible 
choices of Xi and Xi+\. The master equation for the evo- 
lution of probabilities of configurations, is 

dP(C)/dt = -J2TccP(C) + J2 t cc>P(C) . (1) 

c> c 

Since there are X distinct choices for the C into C and 
also for the C out of C, each having probability 1 jX, the 
probability distribution P( Co = a) = cons, independent 
of a, is invariant in time. Thus the probability distribu- 
tion of states on the attractor is a product measure, with 
each site independently occupied with one particle with 
probability 1/2, otherwise being empty. 

In a recent work, Dhar [^9) has shown that the stochas- 
tic Manna model also exhibits the Abelian property and 
is a special case of the Abelian Distributed Processors 
Model. Correspondingly some of the the analytic tech- 
niques of the BTW model also apply to the stochastic 
Manna model. It is only necessary to realize that for 
the stochastic models, instead of associating probabili- 
ties with each toppling, we can assign to each site an 
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infinite stack of random numbers, uniformly distributed 
between zero and one, say. The quenched random num- 
bers in each site's stack then determine the allocation of 
grains during each toppling event. Thus, the q'th ran- 
dom number at (x,t) determines at the q'th toppling of 
that site where the grains will go. There is a one-to-one 
correspondence between any realization of the dynamics 
of the stochastic model, and the dynamics of a determin- 
istic system with a random array (chosen appropriately 
to model the probability distribution of grain allocation) , 
under the same condition of particle additions. 

If we specify the height configuration of the sandpile as 
well as the infinite stack of random numbers at each site, 
Dhar shows that the model is also Abelian. It is easy 
to check that given any unstable configuration with two 
or more unstable sites, we get the same configuration by 
toppling at an unstable site i, and then at unstable site 
i' , as we would get if we first toppled at i' and then at 
i, if the same list of random numbers in the array is pro- 
vided. Iterating this until a metastable state is reached 
proves the Abelian property of the model. 



B. New Results 

The directed stochastic model is also equivalent to a 
deterministic directed model with an infinite stack of 
quenched random numbers at each site. Since the lat- 
ter model is Abelian we can choose to relax each row, 
one site at a time, until it is completely stable, before go- 
ing on to the next higher row. In this case, it is easy to 
see that the deterministic model with quenched random 
numbers shares the same property of invertibility as the 
DR model @. 

Let Co be a starting configuration and R(x, t, q) be the 
infinite array of random numbers, with the initial point- 
ers qoix, t) = for all entries [x, t). The pointers q in the 
array R will move as the sequence of topplings proceeds. 
The i'th particle being added at site Xi and the current 
pointers qi-\{x,t) in the array R known, this results in 
a new stable configuration Cj , and a new set of pointers 
qi(x,t) in the fixed array R. Invertibility follows. In this 
case we are given the current configuration Ci, the cur- 
rent set of pointers Qi(x, i) in the fixed array R and Xi. 
In order to prove invertibility we must determine both 
d-i and Qi-i(x, t). 

On the top row C$_i differs from Ci only at the site 
Xi, with z(xi, 0) in Cj_i being less than its value in Ci by 
l(mod2). If z(xi,0) = 1 in Ci, then no toppling occurred 
and Ci-i is the same as Ci at all other sites; also the set 
of pointers {qt-i = qi}. If z(xi, 0) = then one toppling 
occurred at that site. We locate the pointer qi(xi,0) and 
move it back one step in the stack R{xi, 0, qi(xi, 0)) giv- 
ing qi_i(xi, 0) = qi(xi, 0) — 1. This pointer now points to 
a number that tells us where the two grains were placed. 



The heights at the sites in the second row t = 1 in con- 
figuration Ci-i are the same as those in Ci except at 
the forward neighbors from x\ that received a grain ac- 
cording to R(xi, 0, qi-i(xi, 0)). If both sites received a 
grain then we apply the same procedure to those sites 
as we applied to (xi,0). If one site receives two grains 
then that site must have toppled once. Its height in the 
previous configuration is the same as its height in the 
current one, and its pointer is moved back by one unit, 
determining which downstream neighbors receive grains. 
One continues in this fashion increasing the row t. 

Unlike the DR model, eventually one can encounter 
a site receiving three or more grains from sites in the 
previous row. If the total number of grains received 
at a site, n, is even then the site must topple exactly 
n/2 times. The pointer at that site is moved back n/2 
steps, so = qi(x,t) — n/2, reading the inter- 

vening numbers in the stack at that site to determine 
where the grains from that site are sent. If n is odd and 
in Ci the height is one, then the site must have toppled 
(n — l)/2 times, with its height in Cj_i being 0. Thus 
qi-i(x,t) = qi(x,t) — (n — l)/2. Similarly if n is odd and 
in Ci the height is zero, then the site must have toppled 
(n + l)/2 times, with its height in Cj_i being 1. Thus 
qi-i(x,t) = qi(x, t) — (n+l)/2. One reads the intervening 
sequence in the array of random numbers for that site to 
determine how many grains each downstream neighbor 
receives, and so forth. Thus, given Cj, X4, and qi{x,t), 
with a fixed array R{x,t,q), we can uniquely determine 
Ci—i and qi-i(x, t). This proves the invertibility of the 
dynamics of the SDM. 

For a given array R and set of pointers {qi}, for any 
state Ci there are precisely X distinct choices of Q_i 
and Cj+i corresponding to the X possible choices of Xi 
and Xj+i. It then follows, as before, from the master 
equation for the evolution of probabilities of configura- 
tions that the state prepared with a uniform distribution 
over all stable states is invariant in time. 

Thus for the directed Manna model, the self-organized 
critical state is the set of all stable states with equal likeli- 
hood; it is a product measure state, where the probability 
for a site to be empty is equal to the probability for it 
to have one grain, which are both equal to 1/2. This is 
exactly the same as in the DR model, so for the SDM 
the presence of multiple topplings does not lead to any 
correlations in the states on the attractor. 

IV. DISCRETE EQUATION FOR AVALANCHES 
IN THE CRITICAL STATE 

The fact that the critical state is a product measure 
state leads to a significant simplification; namely the crit- 
ical dynamics can be described as a type of generalized 
branching process. Thus one can simulate or describe 
avalanches in an infinite system as follows. Consider a 
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site which we we will call the origin. The origin in the 
equivalent branching process represents the site that re- 
ceives a grain in the critical state of the SDM. The height 
at that site is either one or zero with equal probability. 
Add one grain to it. If the height now is greater than 
one it topples. Then define the heights at sites (1,1) and 
(-1,1); they are one or zero with equal probability. They 
receive grains from the origin according to the stochastic 
rules of toppling in the directed model, and topple if they 
are unstable. In this way, one can always construct the 
lattice as the avalanche propagates and one can simulate 
the infinite system, albeit always for a finite time. 

We define the quantity n{x, t) to be the number of 
grains added to (a;, t) given that one grain was added to 
the origin. The total number of grains that leave a site 
n ou t(x,t) can at most differ by one from the number of 
grains going in. If n(x,t) is even then n out (x, t) — n(x,t). 
If n(x, t) is odd, then, since the number of grains which 
can leave any site is always even n out (x, t) = n(x, t) ± 1. 
The process is critical and the increase or decrease occur 
with equal probability. Thus we observe there is a source 
of bounded, nonconservative noise, with a threshhold, in 
the dynamics of n during an avalanche that comes from 
the presence or absence of grains in the metastable states. 
Since the number of grains going into a site can only 
arise as a consequence of grains going into its immediate 
upstream neighbors we arrive at the following discrete 
equation 

n{x, t + 1) = ^ (n{x - 1, t) + n(x + 1, i)) 

+0 o (n(x + 1, t))n(x + l,t) + 6 (n(x - 1, t))n(x - 1, t) 
--j{x + l,t)+j{x-l,t) . 

(2) 

On average each site will get 1/2 of the grains going 
into its upstream neighbors. There are two sources of 
stochastic variations from the average. One is conserva- 
tive: Each upstream neighbor may divide its out flowing 
grains unevenly between its two downstream sites, but 
what is taken away from one downstream neighbor is 
added to the other according to the binomial distribu- 
tion. This gives a stochastic current j which is either 
directed to the right (here defined as positive) or to the 
left (here defined as negative) for each site. The first two 
moments of the stochastic current of flowing grains are, 
from the binomial distribution, 

(j(x,t))=0 

(j{x, t)j(x', t')) = n(x, t)5 (x, x')S (t, t') . (3) 

Since this is a discrete equation, the Kronecker delta 
functions, <5 Q , are defined on the set of integers. 

The nonconservative noise is the most interesting and, 
as we shall see, relevant noise. It is associated with the 
fact that the metastable states either add or absorb flow- 
ing grains from the avalanche. However, as mentioned 



before, the number of flowing grains can only change by 
one unit irrespective of the local number of flowing grains 
as long as it is nonzero. This gives rise to the discrete 
Heaviside step functions in Eq. 2 defined as o (u) = 1 
for u = 1,2,3, ••• and 8{u) = otherwise. With this 
convention the nonconservative noise is at each point in 
space-time either ±1 with equal probability 1/4 or with 
probability 1/2. Thus 

(n(x,t)) = o 

(i 1 (x,t)n(x',t')) = ±5 (x,x>)6 (t,t') . (4) 

The appropriate initial condition to describe the 
avalanche is n(x,0) — S o (x,0). The avalanche propa- 
gates and spreads out; eventually it dies out. Then a 
new avalanche, represented by a new realization of the 
branching process is started. 



V. CONTINUUM EQUATION FOR THE 
AVALANCHES 

One could consider a rigorous derivation of the contin- 
uum limit of Eqs. 2-4, taking the lattice size in space, 
A x , and time A t , as well as the grain size, A„, to zero. 
Instead, here we invoke the usual "hand-waving" , coarse 
graining procedure to obtain a smooth function n of con- 
tinuous variables x and t. Expanding to leading order in 
gradients, and time derivatives, we arrive at 

= \ V 2 n(x, t) - 2^M + 29(n(x, t))n(x, t) , 

(5) 

where the threshold function 0(u) = for u < and 
9{u) = 1 for u > 0. By the central limit theorem, the 
noise terms are both Gaussian with first and second mo- 
ments 

(V(x,t))=0 
(n(x,t)r,{x' ,t')) = ^S(x-x')S(t-t f ) 
(j(x,t))=0 

(j(x, t)j(x', t')) = n(x, t) S(x - x')S(t - t') . (6) 

The appropriate initial condition for the avalanche is 
n(x,0) = 5(x). The avalanche grows by increasing or 
decreasing n locally where n is non-zero. Eventually the 
avalanche dies and n(x, t) — everywhere. This equation 
describes the transient out of an absorbing state associ- 
ated with the avalanche. In particular, the state with no 
flowing grains n(x,t) — for all (x,t) is stable, which is 
a requirement of any equation describing avalanche dy- 
namics. 
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A. Analysis 



1. The Threshold Term 



Dimensional Analysis is the simplest tool we can apply, 
and the first step in any theoretical analysis. The dimen- 
sion of the conservative noise is [j] 2 /[a;] 2 = [x] 3 ), 
and dimension of the nonconservative noise is [rj\ 2 — 
(\/[x][t\). Thus, as long as [n] < [x] 2 then the conser- 
vative noise is irrelevant with respect to the nonconser- 
vative noise. Ignoring this term we arrive at our main 
result 



dn(x, t) 
dt 



1 



V 2 n(cc, t) + 20{n(x, t))r){x, t) 



In the region covered by the avalanche n(x,t) > 0, the 
threshold function, 9 — 1 and may be ignored, result- 
ing in an avalanche dynamics described by the linear 
Edwards- Wilkinson equation (27J. Dimensional analysis 
then gives the correct scaling of various quantities. Thus 

1 1 z 

for the SDM, x c ~ t c ' with z — 2 precisely as in the 
DR model. However, the Edwards- Wilkinson equation 
gives a rough surface in one dimension and the maximum 
number of topplings scales as the transverse extent of the 
avalanche . This differs markedly from the 

DR model where n c = 1, independent of the transverse 
extent, x c . 

Continuing with our scaling analysis, the area covered 
by the avalanche is a ~ x c t c ~ t 3 J 2 (as in the DR model), 
but the size of the avalanche includes the extra effect of 
multiple topplings. The size scales as S TLX^tc r ^ tc 
Since on average for every grain added one grain must be 
transported the entire length of the system to the open 
boundary, we have that < s >= T. Since all the geomet- 
ric quantities associated with avalanches exhibit scaling, 
it is reasonable to assume, and can perhaps be proven, 
that the distribution of avalanche sizes, times, and spatial 
extent are power laws, namely P(s,T) ~ s~ T f(s/T D ), 
P t {t,T) ~ t- Tt g{t/T), and P x {x,T) ~ X - T *w(x /T l < lz ). 
The constraint on the average size then gives 1 = D(2— r) 
or r = 10/7. Similarly, from conservation of probability, 
T t - 1 = D(t - 1) = (t x - l)/z, gives T t = D = 7/4 and 
t x = 5/2. 

It is straightforward to check that Eqs. 2-4 also apply 
to the case where there are d transverse dimensions. The 
only factors that are changed are various constants. Ap- 
plying dimensional analysis, we find that in d transverse 
dimensions, x c ~ tc , so that z — 2 and n c ~ x~z~ . 
The upper critical dimension is d c = 2 above which the 
maximum number of topplings does not diverge with the 
size of the avalanche, and mean field results obtain. This 
corresponds to the fact that the surface described by the 
EW equation is flat above two dimensions rather than 
being rough. For d < 2, s ~ n c xft ~ f3/2+d/4^ j e 
D = 3/2 + dj '4 and the other exponents are obtained 
via the above scaling relations giving r = 2 — (4/ (6 + e?) ) , 
and n = 3/2 + d/4. 



Outside the region covered by the avalanche, the 
threshhold function 8 has a major effect on the dynamics. 
In particular, in regions where n{x, t) = 0, the interface 
is pinned and cannot move. The noise does not act where 
there are no flowing grains! This is completely different 
than the usual models of stochastic interfacial growth. 
The threshhold function importantly breaks the transla- 
tional symmetry of the EW equation (n — > n + const.) 
and leads to an absorbing state. Typically absorbing 
state phase transitions have been considered where the 
amplitude of the noise depends on the activity n to some 



positive power |31|. Here we find a very weak effect sim- 
ply distinguishing between having activity and not hav- 
ing it in terms of a threshold function. This effect is 
so weak that the scaling dimensions of the propagating 
avalanche are the same as the linear EW equation. Ob- 
viously if the threshhold function 9(n) were replaced by 
n a in Eq. 5 that would no longer be the case. Thus 
Eqs. 5-6 is a hybrid combining interface dynamics (the 
number of topplings of the avalanche being the interface) , 
and an absorbing state model. It describes a previously 
undiscovered absorbing state phase transition. 



2. Averaging over Noise 



Averaging over avalanches corresponds to averaging 
Eq. 7 over noise and we arrive at a linear diffusion 
equation for the average propagation of flowing grains in 
response to a single grain being added at (0, 0) to the 
critical system: 



d(n(x, t)) 
dt 



1 



whose solution is 



(n(x,t)) 



V 2 (n(M)) 



-x 2 /2t 



(2nt) 1 / 2 



(8) 



(9) 



Obviously, this solution has the important property of 
conservation, namely J dx(n(x,t)) = 1 which is required 
for stationarity. Note that the DR model also obeys ex- 
actly the same equation for the average propagation of 
activity. This equation is enforced by the local conserva- 
tion and symmetry properties of the system and is in no 
way related to the presence or absence of multiple top- 
plings. Thus we get that for all models with the same 
symmetry and conservation of grains, the dynamical ex- 
ponent z = 2 and the exponent r] = 0, since the average 
amount of activity remains constant in time. 



G 



3. The Fokker-Planck Equation 

Ideally one would like to determine the full probabil- 
ity distribution for the number of topplings n(x, t) in 
avalanches. The dynamics of this probability distribu- 
tion P[n;£] is expressed by the Fokker-Planck equation. 
The Fokker-Planck equation can be obtained by straight- 
forward means from the Langevin equation (Eqs. 6,7). 
It is 



dP[n;t] 
dt 



i/„4{,v^p} 



dx- 



Sn 



e(n){^(e(n)P)}\ . (10) 



Unfortunately we are not currently able to analyze this 
equation in any significant way. 



VI. OUTLOOK FOR FUTURE WORK 

A major limitation of the present work is that applies 
only to a set of directed models where all stable states are 
equally likely. Unlike most models of SOC, there are no 
spatial or temporal correlations in the metastable states 
on the attractor. Even in this drastically simplified set- 
ting, the occurrence of multiple topplings has a profound 
affect on the critical properties of the system, changing 
the universality class. This fact suggests that any rea- 
sonable theory of avalanche dynamics in sandpile models 
of SOC must treat the effect of multiple topplings. 

Avalanches are described by the dynamics of particles 
which exhibit an absorbing state phase transition. As 
shown in our main result, Eq. (7), this transition has the 
new feature that the amplitude of the noise is a thresh- 
hold function of the local activity rather than a being a 
power of the activity. This threshhold form of the noise 
amplitude has never before been discussed in the litera- 
ture. It may be interesting to examine a broad range of 
absorbing state phase transitions with such a threshhold 
noise amplitude, as well as the implications of the thresh- 
hold noise amplitude on avalanche dynamics in SOC. 

The picture of avalanches as reaction-diffusion systems 
with an absorbing state was first suggested in ]3^| as ap- 
plicable to SOC systems and later in J33]. In the case 
discussed here, the particles, representing topplings, are 
known to propagate in an uncorrelated environment be- 
cause the probability distribution of metastable states on 
the attractor is described by a product measure. In the 
general case, there will be important correlations from 
the background that must be included along with bound- 
ary effects, leading among other things, to correlations 
between avalanches. It seems to us that Dhar's develop- 
ment of an operator algebra for stochastic models might 
provide a fruitful avenue to pursue further research. 



Note Added: After submitting this manuscript for pub- 
lication, another preprint appeared J34j. They also ob- 
tain the critical exponents for the stochastic directed 
model we discuss. 

We thank D. Dhar for helpful comments on the 
manuscript. This work was supported in part by NSF 
grant DMR-0074613. 
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